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O ' Abstract 

in 

I We extend the quasilocal (metric-based) Hamiltonian formulation of general rel- 

ativity so that it may be used to study regions of spacetime with null boundaries. 
In particular we use this generalized Brown- York formalism to study the physics of 
isolated horizons. We show that the first law of isolated horizon mechanics follows 
directly from the first variation of the Hamiltonian. This variation is not restricted 
bX). to the phase space of solutions to the equations of motion but is instead through the 

space of all (off-shell) spacetimes that contain isolated horizons. We find two-surface 
integrals evaluated on the horizons that are consistent with the Hamiltonian and which 
define the energy and angular momentum of these objects. These are closely related 
5^ I to the corresponding Komar integrals and for Kerr-Newman spacetime are equal to 

the corresponding ADM/Bondi quantities. Thus, the energy of an isolated horizon 
calculated by this method is in agreement with that recently calculated by Ashtekar 
and collaborators but not the same as the corresponding quasilocal energy defined by 
Brown and York. Isolated horizon mechanics and Brown- York thermodynamics are 
compared. 



1 Introduction 

In the almost 30 years since its initial discovery by Bekenstein |l|] and Hawking p[, black 
hole thermodynamics has been one of the most active areas of gravitational research. There 
have been several reasons for this interest, but perhaps the chief one is that these laws 
represent one of the few solid links between the classical world of general relativity and 
the much-sought-after theory of quantum gravity. Just as quantum mechanics gives rise to 
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the statistical mechanics that in turn explains the thermodynamics of condensed matter, 
it is believed that the laws of quantum gravity will explain the statistical mechanics and 
therefore thermodynamics of black holes. The advances made by string theory and canonical 
quantum gravity towards this goal in the last few years are well known (see for a review 
and references). 

That said, there are several ways in which standard black hole thermodynamics differs 
from the more familiar thermodynamics of everyday materials. Perhaps the most obvious 
problem with black holes comes from their very definition. While one can easily localize 
standard thermodynamic systems such as steam engines in spacetime, this is not the case 
for a black hole. By the standard definition, a black hole is a region of spacetime from 
which no null or timelike signal can ever escape. Thus, one must know the entire history of 
a spacetime before one can say with certainty whether or not a particular region is inside 
a black hole. Such a definition is too unwieldy to use except in the very special case of 
stationary spacetimes and so one would like to find a more general class of objects that 
retain the properities of black holes that are essential for thermodynamics, while at the 
same time being more convenient to identify and study. 

A second problem with black hole thermodynamics arises in the definition of thermody- 
namic quantities such as energy or angular momentum. While the entropy of a black hole is 
measured at the event horizon itself, in the usual formulation of black hole thermodynamics 
energy and angular momentum are defined by ADM measurements taken back at infinity. 
Further, to define the temperature of a black hole one must make reference to Killing vector 
fields at infinity. Again this is a situation that is alien to less esoteric thermodynamics. One 
can define the energy of a standard thermodynamic system such as an ice-cube without first 
having to retreat to the outer reaches of the universe. A local definition of energy is obvi- 
ously important too if one doesn't want to include stray galaxies sitting megaparsecs away 
from a black hole in its total mass. Thus, to give a good formulation of gravitational ther- 
modynamics, we need quasilocal ways to both identify black holes (or some generalization 
of them) and to measure all of their thermodynamic properties. 

Quite recently Ashtekar, Beetle, and Fairhurst proposed a solution to these problems 
when they defined isolated horizons Roughly these are closed, non-expanding, null 

surfaces which are in equilibrium with their surrounding spacetime (though as long as the rest 
of the spacetime doesn't disturb the equilibrium of the horizon it may itself be dynamic). In a 
series of papers following the initial two, those authors and others have studied the geometric 
properties of various species of isolated horizons and also shown how they naturally fit into 
the tetrad and spinor Hamiltonian formulations of general relativity. They have provided 
quasilocal, horizon-based definitions of the thermodynamic properties of isolated horizons 
and used the Hamiltonian formalism to derive the zeroth and first laws of thermodynamics 
(see for example for a review of this work and further references). 

Of course, Hamiltonian derived definitions of quasilocal energy and even quasilocal ther- 
modynamics are nothing new and there has been an active literature on these subjects for 
some years (see for example |^, Q and the references contained therein). In fact, the quasilo- 
cal energy community has generally been more ambitious than the isolated horizon one and 
sought to define the energy of arbitrary regions of spacetime rather than just black holes. 
Unfortunately, the various definitions tend to disagree on their localizations and no consensus 
on what is the "correct" quasilocal energy (or even if such a thing exists) has been reached. 
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Nevertheless, there are well developed quasilocal formalisms out there, and it is of interest 
to see what they have to say about isolated horizons and how their quasilocal definitions 
relate to the isolated horizon ones. 

One of the most popular definitions of quasilocal energy and formulations of quasilocal 
thermodynamics was proposed by Brown and York in 1993. Starting with the standard 
Einstein-Hilbert action (with boundary terms) they calculated the first variation of the 
action. Setting it equal to zero they obtained the equations of motion (in the usual way) 
plus a set of boundary conditions. Enforcing these boundary conditions causes the variational 
boundary term to vanish and so ensure that the variational principle is well-defined. The 
conditions vary according to the original boundary terms of the action and so there is quite 
a bit of freedom in the whole procedure. That said, for each choice of the action, one can 
apply a Legendre transform, obtain a Hamiltonian defined on spatial three- surfaces, and so 
define the quasilocal energy (QLE) of any finite section of that hypersurface as the numerical 
value of the Hamiltonian evaluated over the region. The trick comes in picking the correct 
action and it turns out that even in the most simple cases, there is some ambiguity in the 
definition of the QLE. This ambiguity is interpreted as the freedom to pick the zero-point 
of the energy. Similar comments may be made about angular momentum, though there the 
zero-point ambiguity doesn't arise. 

This formalism has many nice characteristics. For example, it reduces to the ADM 
and Bondi [0 energies in the appropriate infinite limits, while in the limit of a sufficiently 
small quasilocal region with volume V, it reduces to the matter energy {V x Tapu'^u^) where 
Tap and u"" are respectively local values of the stress-energy tensor and the timelike normal 
to the spacelike hypersurface |TT| . With the help of Euclidean path integral gravity one can 



use it to (partially) explain black hole thermodynamics (appendix ^ or for more details). 
It has even been used to study gravitational tidal heating However, the formalism as it 
stands is not quite suitable for studying isolated horizons. Specifically, quasilocal quantities 
such as energy and angular momentum are defined by sets of observers who evolve in a 
timelike way. For an isolated horizon, we want to make measurements on the null surface 
itself which means that such observers must have null (or even spacelike) velocities^ 

In this paper, we extend the metric-based Hamiltonian perspective to include null bound- 
ary surfaces (and classes of observers moving with null or spacelike velocities on those 
surfaces). The extension is fairly straightforward from a geometrical point of view, but 
nevertheless the final results differ quite dramatically from those for timelike boundaries. 
In the standard Brown- York work quasilocal energy (sometimes referred to as the canonical 
quasilocal energy or CQLE) is essentially measured by the extrinsic curvature of two-surfaces 
surrounding the system. By contrast, in the new null formulation the QLE of an isolated 
horizon (which we'll denote the null QLE or NQLE) is measured not by the expansion of the 
congruence of null curves that make up the surface (which would be the naive analogue of 
the extrinsic curvature) but rather by the acceleration of those same curves. Further, while 
the CQLE suffers from the ambiguity in the aforementioned zero-point terms, for the NQLE 
this freedom of choice is largely removed by Hamiltonian formalism. This is in agreement 



^The reader might object that physical observers cannot travel with such velocities. That is a valid point. 
Nevertheless from an intuitive point of view it is convenient to continue to think of the boundaries as being 
made up of observers. 
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with the isolated horizon hterature (see for example 0). Finally, the NQLE measured on 
an isolated horizon is NOT equal to the CQLE measured on a timelike surface positioned 
"close" to the horizon. Specifically, in the case of Schwarzschild, the NQLE is m (in agree- 
ment with [^]) while the CQLE is 2m. In fact, for any stationary, axisymmetric isolated 
horizon in asymptotically flat space, the NQLE is equal to the ADM/Bondi mass of the 
entire spacetime. 

In other ways, however, the results correspond much more closely to those of Brown 
and York. In particular, the expression for angular momentum on a timelike surface is very 
closely related to the corresponding expression on a null surface (which in turn we show to 
be eqivalent to the corresponding Komar formula). Further, basically the same Hamiltonian 
boundary terms arise on an isolated horizon as those that show up on a timelike surface, 
though as noted above the roles that they play can be quite different. Also, drawing on the 
extant work on metric-based Hamiltonians, we can, without much difficulty, carry through 
the variational calculations explicitly for the phase space of all spacetime configurations 
rather than just the subspace of solutions to the equations of motion. 

This paper is organized in the following manner. Section ^ establishes notation and sign 
conventions, reviews the definitions and some properties of isolated horizons and considers 
the constants that characterize these objects. The notation will be an amalgam of that used 
in the Brown- York tradition of papers with that used in studies of isolated horizons with 
some new symbols introduced where the systems overlap and clash. Section |^ extends the 
the work on metric-based Hamiltonians to include systems with null boundaries in general 
and rigidly rotating horizons (a close cousin of isolated horizons) in particular. We find 
constraints on possible boundary terms for the Hamiltonian. Section ^, derives the first law 
of isolated horizon mechanics, proposes an integral boundary term for the Hamiltonian that 
both satisfies the constraints and is equivalent to the function found in |TB|, and considers the 



connection between the quasilocal energy /angular momentum and the corresponding Komar 
definitions. Section |] compares the results found for isolated horizons in this paper with the 
Brown- York quasilocal energy and thermodynamics. Section ^ summarizes the results and 
appendix reviews Brown- York thermodynamics. 



2 Set-up 

The calculations in this paper will be done from a Hamiltonian perspective. As such it will 
usually be most appropriate to think not of spacetime, but instead of time- dependent tensor 
fields living and evolving on a spatial three-manifold. At the same time however, many of the 
quantities appearing in the calculations can most easily be understood as four- dimensional 
quantitities associated with three-surfaces embedded in a full spacetime. As such, and even 
though it will create a certain amount of clutter, we will begin by setting things up from 
a three-dimensional perspective and then "stack" the "instants" of time evolution into a 
full four-dimensional spacetime. Only then will we define the various sub-species of isolated 
horizons and then proceed on to the variational calculations. Note that throughout this 
paper, we'll use units where G and c are dimensionless and equal to unity. 
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Figure 1: The region with its boundaries and normals. 
2.1 The 3D perspective 

We begin with the three-dimensional-space plus time perspective, viewing fields as defined 
in a three-dimensional manifold S and dynamic with respect to a time coordinate t. To em- 
phasize their three-dimensional nature we label them with Latin indices. Figure illustrates 
this section. 

Let E be a three-dimensional manifold with Euclidean metric hab, and Ej C S be such 
that Et is topologically S"^ x [0, 1] (where S"^ is the two-sphere and [0, 1] is a real line segment). 
Label the two-dimensional boundaries of this space Aj and Bf and refer to them as the 
"inner" and "outer" boundaries. The induced metrics on these boundaries will be and 
a^fj respectively. Next, if Da is the metric compatible covariant derivative on E then the 
corresponding induced metrics on the boundaries will be and respectively. Often 
where there is no chance of confusion we will drop the superscripts on these quantities. 

Next, define and as the unit normals to the two-surfaces. They are oriented so 
that points back "into" and points "out of" E^. Then, one can define the extrinsic 
curvatures k^b = -|£ncr^ = -cr^^'a^'^Dcn^ and k^^ = -^CnfTab = -(^a^^b'^DcU^ , where 
is used to designate a Lie derivative in the direction v. The reader should take note of the 
sign convention used for extrinsic curvatures. 

The lapse function N and shift vector V"" will define the fiow of time relative to this 
surface in the usual way. That is, Ndt units of proper time pass for every dt units of 
coordinate time, and an observer swept along by the fiow of time will move with velocity 
^JVaV°-/N 0. In particular we will think of the boundaries and Bt as being made up of 
such observers, and so their evolution will be guided by the shift vector field. Thus, given 

^Keep in mind that depending on the spacetime, such "motion" may or may not correspond to "real" 
movement. For example, a particle which is static on a Schwarzschild horizon will have a velocity equal to 
the speed of light, but doesn't actually move by any normal definition of movement. 
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Figure 2: A three-dimensional schematic of M along with its boundaries, normals, and a 
typical foliation surface. No attempt has been made to have the slope of the boundaries 
correspond to their "speed" . 

the position of and Bt at any instant t, the flow of time determines their position in E 
at all future times and therefore the future extent of St as well. 

In addition to the metric field hab-, the Hamiltonian three-surface formulation of general 
relativity requires a surface-momentum tensor density P""^ which is canonically conjugate to 
the metric (and determines the time rate of charge of the metric). Further, for Einstein- 
Maxwell spacetimes, there will also be the vector potential one-form Aq, its canonical conju- 
gate the electric field vector density (which apart from its obvious interpretation measures 
the time rate of change of A^), and the Coulomb potential $ (which joins the lapse and shift 
as a Lagrange multiplier). 

We now embed all of this in a four- dimensional spacetime. 

2.2 The 4D perspective 

Let be a four- dimensional manifold and M (Z M. such that M is topologically [0, 1] x 
[0,1] X S"^. Thus, M is "cylindrical" and bounded by four three-dimensional surfaces Si, 
S2, A, B (see figure Next, let gap be the metric field on M. such that Si and S2 
are spacelike, A is null, and B is timelike. Note that we use Greek indices for these four- 
dimensional quantities. Other situations will be considered later on in the text, but this 
set-up will allow us to capture all of the essential details for now and for later. 

Assume that M is foliated by a set of spacelike surfaces {St : ti < t < such that 
Sf^ = Si and Stj = S2 are the terminal elements of the foliation. These surfaces can then be 
thought of as "instants" in time and in the usual way, we'll say that St^ "happens" before St^ 
if to; < A further orientation is implied by choosing to refer to A as the inner boundary 
and B as the outer boundary. The St then induce corresponding foliations {At : ti < ^2} 
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and {Bt '■ ti <t < on A and B respectively. In the usual ADM/Brown- York formulation 
of general relativity, these foliations would be independent of the spacetime configuration. 
However, that will not be the case here and in section ^ we will see how the Hamiltonian 
formalism restricts the allowed foliations of A (and therefore M) if A is an isolated horizon. 

Using the spacetime metric gap one may define normal vectors, projection operators, and 
metric tensors for the hypersurfaces. Specifically, if is the forward pointing unit timelike 
normal to the then = ~^ u°'up is the projection operator into the surface, while 
= 9aj3 + WqU/j is the iuduced spacelike metric. Similarly, if is the outward pointing 
unit normal to B then = 9a ~ ^afi^ is the projection operator, and 7^/3 = gap — fiafip is 
the induced timelike metric. For the two-surfaces Bt, the projection operator maybe written 
as either of cr^ = h^ — riaTi^ = 7^ + u'^up where n„ is the outward pointing spacehke unit 
normal to Bt in and Ua is the future pointing timelike unit normal to Bt in B. The 
induced metric is then aap = hap — naUp = -^ap + UaUp. Note that 7^ and Ua 7^ Mq 
unless the foliation surfaces are perpendicular to B. That is 77 = u"na = OQ 

Next, let ia be a null normal to A. Thus if A is defined as a level surface of some function 
H of the spacetime coordinates, ia is proportional to daH. Further, assume that it points 
foward in time and into M so that < and i'^na > 0. However, because ia is null, 
it has no natural normalization. To define the projection operator into A, one must invoke 
an auxilarly null normal which is also assumed to be perpendicular to At and which is 
normalized relative to ia so that i'^ia = — 1. Then if Ua is the inward-pointing spacelike unit 
normal to At, one can write 

ia = C.{ua + na) and ia = -^{u" - n") (1) 

for some function ^. That said, the projection operator into A is = g^ + iai^ and the 
induced metric is qap = gap + iaip + iaip- This metric is, of course degenerate and therefore 
qap does not uniquely define its inverse g"^ (though this non-uniqueness doesn't manifest 
itself in any of the quantities that we use in this paper). It should also be kept in mind that 
and ia lie in the tangent and cotangent bundles to A respectively while their compatriots 
and ia do not. The projection operator into At is erf = ga+iai'^ +iai'^ = ga+UaU^ —rian^ ■ 
The induced metric is then aap = gap+iaip+iaip = Qap+UaUp—naUp. Note that aap = qap- 
The metric compatible covariant derivative on M. will be which in turn will induce 
Da as the metric compatible derivative on Sj, on A^, and on Bt- Aq, will be the 
induced covariant derivative on A. Note that it is compatible with qap-, though not uniquely 
determined by it since that metric is degenerate. 

The extrinsic curvatures of these hypersurfaces may be defined as follows. Kap = 
— \Cuhap = —h'^VpU-y is the extrinsic curvature of in M, kap = —crla^pV^ns, is the 
extrinsic curvature of Bt in St, and lap = — ct^ct^V^m^ is the extrinsic curvature of Bt with 
respect to the normal Ua- Each of these last two quantities may be considered to be defined 
for At or Bt as appropriate. Contracted with the appropriate metric tensors these become 
K, k, and I. The addition of a bar over any of these indicates that they are to be calculated 
with respected to the "barred" rather than "unbarred" normal. For example / = —a'^^'VaUp- 

^Note that the sign convention for 77 agrees with |l3 and the transformation law section of [|l5j but is 
the opposite from that originally used by Hayward p6| and the original non-orthogonal boundary paper 
cowritten by the author . 
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On A, the analogous quantity is 9af3 = qZq^^i^s = '^Z'^p'^-y^s, which describes the 
instaneous evolution of the congruences of null geodesies that make up A. The reader 
should take note of the sign convention which is different from that used in the extrinsic 
curvatures. Then, keeping in mind that is normal to A so that 6[ai3] = 0, this may be 
decomposed as 

Oa/3 = 2^Cra/3 + Ca/3 (2) 

where 9 = a°'^6a(3 = -^Cg^/a is the null expansion and (ap = &ai3 — \(^<^ap is the shear. If 
dafB = 0, then Qajj has a unique compatible covariant derivative Aq, on A. 

The flow of time in this manifold is defined by the vector field which is chosen to 
be compatible with the foliation {S^} so that T°'dat = 1. On the hypersurfaces it may 
be written as T" = Nu"' + V"' where again and are the lapse and the shift. It is 
not restricted to being timelike. In fact, near or on A we will expect it to go null or even 
spacelike. Because of the restriction made in the previous section that the evolution of 
the boundaries by guided by the time evolution, we know that T° may be written as an 
element of the tangent bundles TA and TB at the inner and outer boundaries respectively. 
Specifically, on A we have T° = Ni" + where N = N/^ and V" = a'^V'^, while on S, 

= iVw" + where N = A^A, A = /T+V, and (again) V"" = a^V^. In section | we 
will see that the allowed forms of are restricted at isolated horizons. Again this contrasts 
with the usual ADM/ Brown- York Hamiltonian formulations where T" is an independent 
background structure. 

From this perspective the "instants" are a series of hypersurfaces, rather than a single 
three-manifold on which dynamic fields evolve. It will be useful to note that the hypersurface 
momentum P"''' which was an independent variable in the previous section, is closely related 
to the extrinsic curvature Ka/B and can be written as 

The electromagnetic field in Ai is described by the vector potential Aa which defines the 
electromagnetic field tensor J^a/s = da A/3 — df^Aa- Then, the three-dimensional vector po- 
tential is Aa = h^A/3, the Coulomb potential is $ = —AaU", and the electric field density 
is = —\/h/{ATr)J-'"^Ui3. Throughout this paper we will assume that a single Aa can be 
used to cover the entire region Ai. This assumption means that the electromagnetic U{1) 
gauge bundle will be trivial and therefore there are no magnetic charges in - including 
magnetically charged black holes. Results for magnetically charged spacetimes can be found 
using electromagnetic duality. 

It should also be kept in mind that constant t surfaces in Scwharzschild/Kerr spacetime 
in the usual static/stationary coordinates are NOT examples of the surfaces that we are 
considering here. As is well known, such constant t surfaces go null at the event horizon and 
timelike inside of it. Our are restricted to remain spacelike and so one should instead 
think of a foliation such as that induced by Painleve-Gullstand-type coordinates. These 
are analogous to Eddington-Finkelstein coordinates except that they are tied to infalling 
timelike rather than spacelike geodesies. Descriptions of these in spherically symmetric and 
Kerr spacetimes may be found in ||T9[ and |^ respectively. 
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2.3 Review of isolated horizons 



This section reviews the definitions of the isolated horizon sub-species as well as some of 
their properties. It is largely based on the relevant sections of and and the reader 
is directed there for further discussion and derivation of the definitions and results. 



2.3.1 Definitions 

Definition 1 A null surface A is a non-expanding horizon if: i) its null expansion 9 is 
zero everywhere on A for any null normal ii) all the equations of motion hold at A, and 
Hi) the stress tensor Tap of matter fields is such that — T^£" is future directed and causal for 
any null normal . 

A couple of notes are in order on these conditions. First, it is trivial to check that if 
conditions i) and iii) hold for any null normal vector field £° then they will hold for all such 
vector fields and so the definition isn't based on the normalization of Second, condition 
iii) is implied by the dominant energy condition which in particular holds for the Maxwell 
fields that we consider in this paper. Thus, from our point of view, the only real restriction 
in the above is that 6* = 0. Note that any Killing horizon is automatically a non-expanding 
horizon. 

From this simple assumption, many properities follow. First A may be thought of as a 
congruence of null geodesies. By definition they are surface forming and so their twist 9[ap] is 
zero. With the help of the Raychaudri equation one can quickly see that ^ = then implies 
that their shear must be zero too. Therefore Oap = 0. It then follows that the induced 
Ao is the unique compatible covariant derivative for qap. From the same calculation the 
twice contracted stress-energy tensor T^^^la^p = as well. In turn, T^^laC-p = implies that 
fa = a2,F^si^ and /„ = a^^F^si^ are both zero (or equivalently the electric and magnetic fields 
projected into Aj are zero). Physically this means that there is no fiux of electromagnetic 
radiation across the horizon (the relevant component of the Poynting vector is zero). 

Moving on, direct calculation and Oa/s = implies that there exists a one form Ua such 
that qZ'^-yi'^ = ujai^- Explicitly, 

uJa = -f^eia + (y^i^^fil', (4) 

where kh = —iai^^iB^'^ is the acceleration of the null normal Note that Ua is defined 
entirely with respect to quantities intrinsic to A. These quantities are, however, dependent 
on the normalization of and under a redefinition of the null normal = fi", 

Ki: = fKi + Cif and iL)'^ = Ua + qi^pf- 

Definition 2 A non-expanding horizon A becomes a weakly isolated horizon (A, [£]) if 
it is equipped with an equivalence class of normals [£] such that q^CeUp = for all i G [i] 
(where i' ^ i if i' = ci for some constant c) and a gauge choice such that q^CiAp = 0. 

It quite easy to show that for any non-expanding horizon A there are many choices of [i] 
and gauge choices for the U{1) gauge so that these requirements are satisfied. Further, given 
such choices, ki and $f = —i'^Aa are constants over A. However, even with the equivalence 
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class chosen it is easy to see from the transformation law given above that the value of K£ 
isn't fixed. A specific G [i] must be chosen to fix that constant over A. Similarly, enough 
gauge freedom remains so that may be assigned any value that is desired (though once 
chosen that value is fixed everywhere on A). 

Definition 3 A rigidly rotating horizon (RRH) (A, [i], if), is a weakly isolated horizon 
with rotational symmetry. That is, L^pi'^ = 0, CipQais = 0, C^^Ua = 0, C^i^q^Ap) = 0, and 
C^{q2,qpF.ys) = 0. Further (p" has closed, circular orbits, and is normalized so that those 
orbits have affine length 2tt. 

We will say that a specific foliation A^ is adapted to an RRH if G TA^ for all t surfaces. 
Correspondingly, the time flow field is adapted if T" = — fii^y?" for some G and 
some real constant Vl^^. Then, an adapted flow of time selects a specific "normalization" of 

and therefore an acceleration ki as well. Equivalently, we can say that a choice of and 
VL^p selects a T" on A. In that case T° is no longer an independent vector field but is instead 
defined with respect to quantities intrinsic to the surface^. 

The canonical example of a rigidly rotating horizon is, of course, the event horizon found 
in Kerr-Newman spacetime. Any linear combination of the timelike and angular Killing 
vector fields of the spacetime is then an adapted flow of time. Keep in mind however, that 
the t = constant surfaces in the standard Boyer-Lindquist coordinates are not spacelike 
everywhere and therefore not a suitable foliation from the current perspective. Instead one 
should think in terms of the Painleve-Gullstand-type coordinates defined in [Q. 

Finally, for completeness we note that an isolated horizon is a weakly isolated horizon 
on which and Aq, commute. We won't actually make use of these in this paper, though 
often will follow convention and use "isolated horizon" as a generic term for any of the 
surfaces discussed above. 

2.3.2 Constant quantities on the horizons 

Given a non-expanding horizon A which is foliated by {A^ : ti < t < ^2} one can define an 
well as the electric and magnetic charges of each surface A^: 

A^it) = [ d^x^, (5) 

Qi^it) = f d^x^Ey, and (6) 

J At 

PA(t) = I d^Xy^B^. (7) 
J At 

In the above, E\- = {l/4TT)Fai3i°'i^ = {l/4:TT)Fai3u"n°' is the normal component of the electric 
field to At, while B\- = —{\/ATx)F*pt l^ = — (l/47r)F*^M°n" is the normal component of the 

^This may set off warning signals for readers who are familiar with metric formulations of the action 
principle. In those cases, St and T" are taken as background structures independent of the spacetime 
variables and therefore ST" — 0. Thus one might worry that this change will disrupt the variational 
calculations of the next section. This is not the case as they are based on Hamiltonian rather than Lagrangian 
variations, and so don't assume that ST" = 0. As far as the Hamiltonian is concerned T" doesn't exist, so 
whether it is varied or not is a matter of indifference 
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magnetic field. F*^ is the dual electromagnetic field tensor. Note too that though we have 
chosen to write Qa* and Pa^ in terms of the normals and field tensors, they could also have 
been written purely as integrals of the pull-backs of the field tensors into Aj in which case 
their independence of the normals is manifest. 

Now, it is quite easy to show that if the At are all closed surfaces (which they are 
by construction), then these quantities are constant in coordinate time and therefore not 
dependent on which surface A^ they are evaluated on. Therefore, the area and electric and 
magnetic charges of a non-expanding horizon are invariants of A and independent of the 
choice of foliation (though of course we have already set Pa = with our assumption that 
a single Aa covers 7V1). 

Next, moving up a step to weakly isolated horizons and continuing to assume that the 
At are closed, one can show that for any constants ci and C2 and vector field f " G TA such 
that £iv'' = 0, 

CAiv) = [ d^Xy/av''{ciUJa + C2E^Aa), (8) 
J At 

is a constant in coordinate time, and therefore an invariant of the horizon. In particular, for 
a rigidly rotating horizon we can define the invariant 

Ja = - f c/'xv^<^"(c^„/(87r) - PkA), (9) 

J At 

which in section ^ will be seen to be the angular momentum of the horizon. 

Thus for a rigidly rotating horizon we have four geometrically defined invariants: A a, 
J A, Qa, Pa (though we have a priori assumed that the magnetic charge is zero), ne and fl^ 
are also constants over the surface if we single out a specific adapted flow of time. Later 
we will see that in order for the Hamiltonian formalism to be well defined, these constants 
(along with ^i) will have to be functions of A^, Ja, and Qa (though there will be quite 
a bit of freedom in defining the functions). Therefore, will be a "live" vector field and 



dependent on the spacetime over which it is defined 
Finally, for a given adapted time flow, 

Ea{T) = J^d^x^T''(^ujJ{A7r)~E^{2Aa-^iL)) (10) 

= -^Aa + 2fi^ Ja + $^<5a 
An 

is also an invariant. If E^iT) could be interpreted as the energy of an RRH then this would be 
a Smarr relation for RRHs connecting their energy, angular momentum, and electric charge. 
In fact, in section ^ we will see that Ea{T) is very closely related to the Komar/Bondi energy 
for a stationary spacetime and so can be taken as the the total energy of the RRH. 

With this review under our belts and promise for the future in our minds, we turn 
examine the quasilocal Hamiltonian formulation for general relativity in the presence of null 
boundaries. 
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3 The metric-based Hamiltonian 



Over a four-manifold Ai with no boundary, the standard action for gravity coupled to 
electromagnetism is 

Iem = / d^'x^in - A - (11) 
lovr Jm 

TZ is the four-dimensional Ricci curvature scalar, A is the cosmological constant, and as 
defined earlier, gap and F^p are the metric and Maxwell tensors. Taking the variation of this 
action with respect to the metric and solving 5/ = one recovers the standard equations of 
motion. 

If, however, one considers a manifold M with boundaries, then boundary conditions 
must be enforced on dM in order to keep the first variation well-defined. It is well known 
that these conditions encode themselves in the variational principle in the form of boundary 
terms added to the standard action (see for example |^ for a discussion of this point). These 
boundary terms are well understood for timelike and spacelike boundaries but not so well 
studied if a boundary is null. Thus, an important part of the work of this paper will be 
to derive the appropriate boundary conditions for such null boundaries. This is done by 
calculating the variation of ([TT]) including the effects of the boundaries and then reading off 
which quantities must be fixed to keep the variation well defined. If necessary boundary 
terms may be added to change which quantities are to be fixed. 

However, at the same time, we are interested in working within a Hamiltonian rather 
than Lagrangian framework. Thus, we want to take the Legendre transform of Iem to get 
the corresponding Hamiltonian. The first step to taking that transformation is to rewrite 
the action in terms of quantities defined on the three-dimensional hypersurfaces rather than 
purely geometrical quantities on a four-manifold. To that end, and assuming the full set of 
boundaries defined in section |2.2| (and where for now A is just a null surface rather than an 
isolated horizon of one type or another), we may rewrite equation ( |Tl| ) as 

Iem = - I d^xiP'^fhap) - I d^x{P^V^) - f d^x{P^^) (12) 
JE21 •'A21 J B21 

+ j dt^j^ d'x{P'''^CTKp + £''CTAS\ 

+ j dtl^J^^ d'xiP^Crv^) + J^^ d'xiP^CrV^)^ 

+ j dt^- d^xT^'iLUa/iSn) + E^Aa) + d^x^T'^iu'V an^/{M) + E^aS^ 
— j dtHuk 

where /^^^ = Js^-Jsi, /a^i = /a^-Jai. and /^^^ = Jb.-Jb,- P% = l/(87r)loge and 
P^ = — l/(87r) sinh ^ rj will be seen to be momenta conjugate to y/a on At and Bt, while 
P"^ and S"' are the hypersurface momentum and densitized electric field conjugate to 
and Aa respectively. Further, 

Huk= f d^xiNH + V^Ha + T'^AaQ) (13) 
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where Ti = —\/h/{S'K){Ga/3 + ^Qap ~ 87rTai3)u°'u^ , and Ha = —y/h/{8'n')h2^{Gyi3 + — 
StiT^p)u^ are the standard energy and momentum constraint equations for spacehke hyper- 
surfaces, while Q = —DaS°' is the Gauss constraint equation. 

From a Hamiltonian perspective each of these quantities may be recast as defined entirely 
on Sj. Then, with some prescience of how the calculations will come out, and which boundary 
terms we will want to fix, we define the new action, 

/ = Iem + [ d?x{P^^Kp) + / Sx{P%^) + / Sx{P^^) (14) 
JS12 •'A21 J B21 

+ [ d^xy/aT^iua/iSn) + E^Aa) - [ d^x^^Q, 

J A JBt 

where B = — 7"^Van/3 is the extrinsic curvature of B in M. With the exception of the 
term at A, these are the standard terms added to a quasilocal action to ensure that its 
first variation is well-defined when the boundary metrics (and potentials) are held constant 



||T6| , [T5|. They are independent of the time flow, and so the action with these boundary 
terms remains similarly independent. By contrast, the term on A explicitly incorporates 
the time flow. As we shall soon see, this inclusion means that the Hamiltonian formalism 
places restrictions on the allowed forms of T". This is in marked contrast to the traditional 
formulation where it is an independent background structure. That said, the new action 
may be rewritten as a coordinate-time-integrated functional of quantities defined entirely in 
St: 

/ = Jdt{J^ d'xiP'^'hab + S'^Aa) + d^x{P%^) + d2^(P^v^)| (15) 

- j dt^^Huk + d\./^{N{e + e'") - V\ja+Ja))] . 

We are now adopting a Hamiltonian point of view, and so the interpretation of many of 
the quantities in the above expression has changed. First and foremost, the momenta P'^^, 
E"", P^, and P^ are now independent fields with no a priori connection to their conjugate 
configuration variables. In particular, the time derivatives of those configuration variables 
(indicated by dots) are no longer geometrically defined by Lie derivatives in the T" direc- 
tion. From the three-dimensional perspective this direction no longer exists and the time 
derivatives are now freely determined as well (though of course the variational principal will 
ultimately link them back to the conjugate momenta). The boundary terms on Bf are well- 



known from the work of Brown- York [Q and its non-orthogonal generalizations ||l^, |15[ . In 
three-dimensional language: 

e = k/{87i\) + 2riP''''nanb/Vh, (16) 

5- = -(^Wv^)('^'/A-^AX), (17) 

ja = ~2aacP'''nd/Vh-\/i87i)daV, ^nd (18) 

: -{S'n,/Vh)alA,. (19) 



-m 
J a 



Their meaning in more clear from a four-dimensional perspective. There, e = k/ (Sir) where k 
is the extrinsic curvature of Bt with respect to the normal n" , = E\-^ where $ = —u^Aa 
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is the Coulomb potential with respect to m", Jq = l/(87r) erf ■u'^V/jn^, and = E^^a^Ap. 
Physically, the e terms are related to energy while the % terms are connected with angular 
momentum [0. 

Applying the Legendre transform, we identify the terms on the first line of ([l5|) as kinetic- 
energy-type terms and the last line as the negative of the Hamiltonian H (in standard 
classical mechanics the analogous relationship is / = — H) . The Hamiltonian equations 
of motion and the corresponding boundary conditions are found by solving 5/ = where 
the action is in its dynamical three-surface form given in equation (0). The only really 
troublesome part of the calculation is finding the variation of H^k. Luckily, however, its 



variation has previously been calculated in |]14| and for the pure gravitational case (using 
slightly different notation) in |18|. Then, varying with respect to hab-. A,, S"', N, and 
V"" we have 

SHuk = I d?x{{n-<^Q)6N + {7ia~AaQ)5V''-NQ6^)] (20) 
+ / <exUN{^ + vn''Aa)6[{V^/^/h)£^na] - ^^[{^ / Vh)£''nbAa]\ 

JBt-At ^ ' 

d^x {N{V^/\^){r + vn,e^''''Bb)SAa} , 



where [hablr, [P"'^]t, [AoIt, and [S"']t are each equal to the corresponding Lie derivatives 
in the direction (though of course they have, in this case, been obtained from the vari- 
ational rather than a geometrical calculation), e = A;/(87r) (the energy density on Bt 0), 
= 2/ \/hP"-^nanb (the time rate of change of ^ 0]); and ja = —2aabP^'^nc/\/h (the 
angular momentum density on Bt [0]). v = V^Ha/N (and can be thought of as the "radial 
speed" of the boundary relative to S^), V"" = s""^ = (l/Svr) (fc'^'' - (k - n'^aja'^'') (the 

stress tensor on Bt induced by the gravitational field [0]), Qa = {l/N)DaN (which from 
a four- dimensional perspective is Oq, = u^VpUa - the acceleration of the timelike normal), 
sl-^ft = (l/Svr) (/"''-(/ + u^al)a"^) (the "time" stress tensor [l|]), and al = n^VpUa (the 



acceleration of the spacelike normal along its length), e"'^'^ is the Levi-Cevita tensor on 
and = Vh/{ATT)e^^^DbAc is the densitized magnetic field. 

Now, if — 1 < f < 1 then the evolution of the boundary Q is timelike, the case that 
was considered in great detail in [T^, [T^, |TB[ (where the reader is directed for further 



details on the result that we will quote below for Bt). On the other hand if f = ±1 then the 
evolution is null and we are considering a situation which has not been directly investigated 



until now (though it was considered as a limiting case in [jTO[ and [0). We will set f = 1 
[v = —1 can then be considered simply by reversing the sign of Ua). Thus the boundary Q 
expands/contracts with the speed of light in the direction (though for the non-expanding 
horizons that we consider that actually means it isn't expanding at all). 

Then harking back to the definitions surrounding equation (^, one can write 

0ab = -Ukab + lab) (21) 
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Ua = -cryKbc + S-Kd^P^ (22) 
= ^{^n^D^N~n^n'K^, + ^iP%-CyP^)^ (23) 

where one recognizes that in the Hamihonian context, ^ = exp(87rPy^), Kab = ^{Phab — 
2Pab), and lab = (^a'^t-^cd- With all of this is mind (and with a quick consulation with [|n| for 
the details on B) we can write 6H (recall that this is the negative of the last line of equation 
(0)) as: 

6H = [ d^x{{n-^Q)6N+{na-AaQ)6V''-NQ6(!?)} (24) 

+ / d'xllKtjTSP"' - [P''']T5Kb + [Aa]T5£'' - [SVAa} 

+ d^x{N^e6[i^/^/Vh)nbS''] - V^/{ATT)Nr6Aa - V''6[i^/^/Vh)n,£'Aa]} 
+ 1^ d^xV^{{e + e^)5N - (j, + J:)5V- - (iV/2)s»''5a„,} 
+ / d^xN^/Vh\\8''na]5^ + Bani,e''^^5AA 

JBt ^ J 

+ / d'x[[^]^5P%-[P%]r5V^} 
+ ld'x[[^]^5P^^-[P^]r5^}, 

J Hi 

where of the newly appearing quantities, s'^* = ■^{0"'^ + [i^e. ~ 6]a'^'^), s"-'^ = -^{k"'^ + [^cn^ — 
k]a-''), = Vh/i^7T)e-''''ua'^bA,, [V^J = -iV^^ + Cy^, [y^^ = -Ny/^T+ Cy^, 
and [-P^]t and [-P^]^ are freely defined functions. 

The Hamiltonian equations of motion and attendant boundary conditions are then be 
found by solving 5/ = 0. We have: 

51 = (Terms that vanish for solutions to the equations of motion) (25) 
+ / d^xiP'''6hab + £''6Aa)+ f d^xP%5^/^+ f d''xP^5^ 

-JdtJ^ d^x |iV5[v^^^/(87r)] + V''6[y^ujJ{8TT)] + ^^^'daab^ 

- JdtJ^ d^x {N<!>eS[{^/^/^/h)nbS^] - v^/(47r)iV/"M„ - V''6[{^/^/^/h)nbS^Aa]} 
-Jdt j^d''x^[{e + e^)5N-{j, + f:)5V'^-{N/2)rHaab} 

- [dt f d^xNy^/y/h{[S''na]6^ + i3anbe"^'=(5ic} , 

J J Bt 

where to tidy up an already formidable expression, we haven't written out the equations of 
motion term. It only goes to zero though if the Einstein-Maxwell constraints 7i = 0, Tia = 0, 
and Q = are satisfied, along with the time-evolution equations P"* = [P'^^Jt and S"" = [S°']t 
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(which encode the Einstein-Maxwell evolution equations), and hat = [hablr, = [-AoIt and 
= [v^T (which establish the connection between the conjugate momenta and the time 
derivatives of the configuration variables). Thus, the full four- dimensional structure and 
equations of motion come from the three-dimensional action and solving 61 = 0. 

These equations of motion are, of course, well understood and nothing new. It is more 
interesting to examine the boundary conditions that must be imposed on solutions so that 
the variation will be well defined. In the first place, by the second line of (^), hab and Aa 
should be fixed on the initial and final surfaces. This fixing of initial and final conditions is, 
of course, standard for a variational principle. 

In the meantime, at each surface Bt, the metric, lapse, shift. Coulomb potential, and 
tangential components of the vector potential should all be prespecified as boundary condi- 
tions. Note however, that fixing these quantities does not fix any thermodynamic quantities 
of interest. As shown by Brown and York, the energy and angular momentum both depend 
on Va^a which is not determined by 70,/?. At the same time, the electromagnetic conditions 
are sufficient only to fix the magnetic charge, which we have assumed to be zero anyway. 
The electric charge is left free. Thus, with these boundary condititions, the variation is well 
defined but the allowed thermodynamic quantities on B are not fixed. Note too that for 
any function/functional Jq of the induced metric 7^/3 and vector potential 7fw4,/3, 6I0 = 0. 
Thus, any such term may be added onto the action without affecting the equations of motion 
or even variational boundary terms. This is the source of the zero-point energy ambiguity 
discussed in the introduction. The obvious choice is to set Jq to zero and then forget about 
it, but unfortunately for such a choice the action is non-zero and even diverges for infinite 
regions of Minkowski space, let alone any curved spacetime. Thus, Iq should be chosen to 
compensate for this problem. A discussion of the various choices that are made may be 



found in [|l^, |14|, but here we just need to keep in mind that the choice is not determined 
by the formalism. 

Now, let us turn our attention to A. If it is an RRH with an (unspecified) adapted fiow 
of time, we can use the properties derived in section |2.3| to show that the boundary term at 
A becomes: 

SIa = -J dt (^5^A + ^iSQa + %6Ja^ . (26) 

Since all of the quantities are constant in time, the integral becomes a multiplicative factor 
of ih-h). 

Now, a simple but at the same time slightly subtle point must be made. For a given 
isolated horizon, we have seen that Aa, J a, and Qa are constants and so at first blush one 
might think that the above term is automatically zero, as for the corresponding the situation 
on B where the fixed boundary metric sets the variational boundary term to zero. On A 
however, we haven't specified that it is a specific RRH, but only that it is rigidly rotating 
horizon. Therefore, the variation considered includes the freedom to move between different 
isolated horizon solutions and consequently different values of those constants. Consequently, 
these terms are not constant with respect to the variation. Thus, as things stand, 61 is not 
zero for variations between these perfectly well-defined solutions to the equations of motion. 
However, this can be fixed as there is still the freedom to add a boundary term on A 
(just as a boundary term has already been added at B in equation ([T5|) in order to fix 
the boundary metric and ensure that no thermodynamic quantities have been inadvertently 
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fixed). Then, Ja must take the form 



Ia = -I dtHAiAA, Qa, Ja), (27) 
where must satisfy the relations 

oAa Stt oJa oQa 

This set of equations only makes sense if k^, and V^^p are themselves functions of A a-, 
Qa, and Ja though the form of their dependence is not fully determined by 61 = 0. As has 
already been seen, the choice of a K£ and a is equivalent to a choice of T" on the boundary 
A. Thus, the allowed are now restricted by the spacetime configuration in contrast to 
the usual quasilocal formulation where the is an independent background structure that 
may be freely chosen. For a discussion of this from the symplectic point of view, the reader 
is referred to [p^, [13] . 



4 Physical quantities and thermodynamics 

We now turn to a physical interpretation of this mathematics. 
4.1 The first law of RRH mechanics 

To make the physical connection with thermodynamics, let us restrict ourselves to the space 
of all (magnetic-charge free) solutions to the Einstein-Maxwell equations that have an inner 
boundary A that is a RRH and an outer boundary B which has a given induced metric jais 
and vector potential 7fw4./3. Further, we will assume that B has a timelike Killing vector 
field and so its physical fields are stationary. As discussed in [|^ this means that the value 
of the Hamiltonian boundary term evaluated on B is time (and therefore exact Bt surface) 
independent. Then, the value of the Hamiltonian for any element of this space is given by 

Hsoi = -Ha{Aa,Qa,Ja)+ f d'xV^(N{e + e"'-e^)-V''{ja + (29) 

though it doesn't matter which surface Bt is used for the evaluation of the outer boundary 
term. We may write the outer boundary term as the functional iJ^fn", ^Q,n"] since the only 
freedom left to it comes from and n^Aa- £^ is an extra boundary term arising from the 
reference term /q. Its exact form doesn't matter here, though keep in mind that it is usually 
chosen so that Hgoi (without the inner boundary terms) matches the ADM energy in the 
appropriate limit for asymptotically flat spacetimes. 

In spite of the above restrictions, the phase space under consideration is still extremely 
large. To wit, the inner boundary term demands that A be an RRH, but does not require 
that it be any particular RRH. The outer boundary term fixes the induced metric and vector 
potential on B, but does not place restrictions on the other components of those quantities. 
Thus even though the outer boundary term defines the ADM mass if A4 is asymptotically 
fiat and B is taken to timelike infinity |^, the variations are such that the value of that 
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mass is not fixed. In the same way, the total electric charge and angular momentum of the 
spacetime are not fixed in that limit. Similarly, as was discussed in the last section, the 
quasilocal quantities are undetermined even if B is finite. Further, though the assumptions 
mean that matter/energy is not allowed to flow across either boundary, there is still a large 
amount of freedom left to such flows inside M. A study of the range of spacetimes containing 
isolated horizons may be found in ||23[ . 



Now given that Hb corresponds to the CQLE/ADM mass inside B it is natural to 
interpret as the energy of the RRH. As it stands though, this is not very satisfactory 
since there is so much freedom left in its definition. This problem will be addressed in the 
next section where we'll investigate the Kerr-Newman region of the phase space of solutions. 
For now though, let's tentatively make this association and see where it leads. 

With respect to variations through the phase space of solutions, we know from equations 
(H) that 

SH^ = ^6 A A + n^SJA + ^iSQa. (30) 
ovr 

Then, interpreting SH as the energy of the RRH, Hi as the surface gravity/temperature, J a 
as the angular momentum, and fl^ as the angular velocity conjugate to Ja, this is the first 
law of thermodynamics. However, note that there is still freedom in the definition of k^, 
fl^p, and the Coulomb potential and so at the moment this is just the first law of isolated 
horizon mechanics and it is important to note that it holds for all of the possible forms of 
K£, Q^, and The next section will make the physical connection but for now we point 
out that this law must hold for the Hamiltonian evolution to be consistent - or conversely 
the Hamiltonian evolution gives rise to the first law. This result was found via symplectic 



methods in [22| and [13| 



It is not hard to show that variation of the outer term in this case is 

6Hb = -Hsn'^dn^, (31) 

and we note that it is not coupled in any way to the variation on the inner boundary. This is 
not surprising given the freedom of the phase space. Energy can be injected into the region 
M without affecting the RRH itself, so there is no reason why the inner and outer boundary 
energies should change in lock-step. An example of such a process would be to insert a 
spherically symmetric shell of matter around a Schwarzschild black hole. Then, the metric 
inside the shell would be not be affected while the ADM energy would certainly increase. 
More complicated and dynamic situations could also be considered. 



4.2 Physical interpretation and calibration 

The quickest way to extract the physical content of the law and the quantities that it connects 
is to consider the Komar expressions for angular momentum and energy. First, given a 
stationary spacetime that contains an RRH and has a global rotational Killing vector field 
Lf"", the Komar angular momentum integral evaluated on a foliation element Aj is 

J A = I <fxV^{ti^-tl^)V^^p (32) 
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which of course is equal to the purely gravitational part of Ja- At the same time it was seen 
in |jl3| that the purely electromagnetic part of Ja is equal to the total angular momentum 
of the attendant electromagnetic fields. Specifically, if B is stationed at infinity, is 
asymptotically flat, and is a global angular Killing vector, then for any surface St, one 
can show that 

J^-Ja=- d'xVhT^^u'^^^ = Jl^' (33) 

J St 

where is the Bondi mass, and Ja*^ electromagnetic component of Ja- Therefore 

Ja = Joo for stationary spacetimes and so measuring Ja at the horizon, one gets not only the 
angular momentum of the hole itself but also the angular momentum of all of its associated 
electromagnetic fields. This is, perhaps, a bit of a surprise since we have tacitly been working 
with the assumption that Ja measures the angular momentum inside At. On the other hand, 
given that the mathematics has viewed A not as a boundary containing a black hole, but 
rather as a boundary of M perhaps this isn't such a surprise afterall. 

Encouraged by this result, let us now consider the Komar integral for the energy. In that 
case, if an adapted time flow vector field T"" = 1°- — fi^y?" can be extended to a global Killing 
vector field, the Komar integral for the energy can be evaluated on the horizon as 

El = ^ \ (fx^i^ti" - fr)VaTf3 (34) 



in J At 

AtT J At 



K 



K 



— . AA + 2f2(pJ^, 

which is equal to the purely gravitational part of the constant E/\ defined back in equa- 
tion (0). Given the previous result for angular momentum, it is tempting to guess that 
the remaining electromagnetic part of -Ea is equal to the stress-energy of the surrounding 
electromagnetic fields. In fact this is the case, for if M is an asymptotically fiat stationary 
axisymmetric spacetime with global Killing vector field T", then one can show that 

E^-El = 2 [ d'xT^pu^T^ = Ei^. (35) 

J St 

That said however, we have not yet shown that E^ is a satisfactory H^. Let us now postulate 
just that, and see where it leads. 

First, if we assume that J^a = E/^ then a Smarr-type formula holds for isolated horizons. 
Thus, let us instead make the slightly more general assumption that 

Ha = + 2%Ja + $,Qa. (36) 

Any consequences of this assumption will then automatically hold for the special case where 
i^A = -E'a as well. Then, equations (^8]) expand to a set of three coupled differential equa- 
tions. Namely 

= / + 7^7^7^ + 27^'^a + 7T7^Qa 37 
Stt a^A 47r aAc^ c/Aa 
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= Tvf 7^ + ^T^^A + ^Qa. 39 

Quite a mess, but luckily these equations decouple when we consider the relationships be- 
tween the partial derivatives of k^, f2<^, and $^ that may be derived from equations (|28|). 
Partial derivatives of i^A in phase space commute so 



1 dKii dQ^ 1 dKi d^i dVL^ d^i 



SndJA OAa' SndQA dA^' dQ^ dJ^' 
and equations (p^),(p8D, and (^) become 



(40) 



= K£ + 2— ^Aa + 2— ^ Ja + t^Qa, (41) 



Ma 5Ja <9Qa 

3^ BQ dQ 

= + 2^Aa + 2^ Ja + ^Qa, and (42) 

C/iA CJa Cv/ 




= 2— ^Aa + 2-/Ja + ^Qa. 43 
oAa oJa oQa 

These have the general solutions: 

t^e = —F=J —F= , = —F=g —r= , and = 

v^4 V^A \/Aa) \fA \Aa -JAa) 

where f{x, y), g{x, y), and h{x, y) are freely defined functions on the real plane. Thus, given a 
Smarr-type formula for Ha (which in particular includes the case where Ha = Ea), then k^, 
Q^, and $^ must have the functional dependences given above for the Hamiltonian evolution 
to be well-defined. In any such case the first law of isolated horizon mechanics will hold. 

Now, that said, for at least one region of the RRH phase space there is a natural choice 
for and therefore ki, Q^, and <I>^. For the Kerr-Newman spacetimes, it is natural to take 

as the stationary time Killing vector normalized to have length —1 at infinity. Then for 
this choice of and also gauging Aa so that $ = at infinity, 

2RiJiRl + Qir + AJi 



% = ; ^'^^ (46) 



^{Rl + Qlf + 4Ji 



where Ra is defined by A a = 47ri?^. It is easy to check that these take the forms required 
by (^4|). Therefore by the Smarr formula 



2R-A • ^^^^ 
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which of course is equal to the ADM energy in Kerr-Newman space (or the Komar mass 
evaluated at infinity). Because they meet the conditions ( ^41) one is perfectly free to extend 
these definitions for ki, VL^p, and across all of phase space and so effectively calibrate the 
invariants so that they will match their natural values in the Kerr-Newman section of the 
phase space. 

This calibration was first made in ||13| from a slightly different perspective which it is 
illuminating to consider. There, no postulate was made about the form of i^A? but instead 
it was simply required that k^, and take the above values across the phase space so 
that they would match their natural definitions in the Kerr-Newman sector. Then, the first 
law was integrated and found to give rise to equation (|^). Thus, from this perspective if 
one calibrates k^, and using the Kerr-Newman sector, then the given form of H^^ 
follows from the first law with no further assumptions. It is a derived quantity and though 
the f2<^, and have been specifically chosen to match Kerr-Newman values, has not 
but rather the formalism has forced it to take that form. 



5 Comparison with Brown- York 

The procedure used to calculate the variations in this paper is a straightforward extension 
of the calculations found in the metric-based Hamiltonian literature. In particular, it is 
very closely related to both the author's PhD thesis [0 and the recent paper by Brown, 



York, and Lau |T8| which depart from earlier works by calculating the variations directly 
from the three-surface based Hamiltonian rather than from the four-dimensional action and 
so obtain Hamiltonian rather than Lagrangian equations of motion (though the boundary 
terms remain the same). Thus, it is of interest to compare the current work with those earlier 
ones in the Brown- York tradition that derive quasilocal energies and thermodynamics for 
regions with timelike boundaries. First, we examine the quasilocal energies. 



5.1 Quasilocal energy 

Let us compare the expressions for the quasilocal energies. On the RRH, the NQLE can be 
written in integral form as 

Ha = <fx^[K,/{ATx) + V^"d;„/(47r) + Ek($, + 21>"A)} (49) 

while for a general Bt in B, the corresponding Hamiltonian boundary term may be written 

as 

Hb = d^xV^ {Ne + Vj^ + E^{2N^ + 2V'' A^) } , (50) 

where for now we have neglected the reference term arising from Iq. Of course for if a, T°' 
is required to be tangent to the null surface A while for Hb, is required to be tangent 
to the timelike surface B, so the timefiow vectors for the two are not the same. That said, 
we can compare them by considering a surface = Bt formed when a null surface A and 
a timelike one B intersect. Then there is a different timeflow vector for each surface, but 
nevertheless we now have the two QLEs defined on the same surface and so can compare 
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their functional forms. To faciliate this comparison, we assume that the each of the timeflow 
vectors projects to the same shift vector V^" in = Bt. . 

We consider the similarities, starting with the V"' dependent terms. First note that the 
electromagnetic angular momentum terms are identical. The gravitational angular momen- 
tum terms are also very closely related. A straightforward calculation shows that 



-aiPVpi^ = a^u^Vpn^ + d^i = 8nj^ + log^ (51) 



Thus, for cases where Cy^ = 0, the angular momentum terms are identical. Of course, given 
the near identical derivations of the two QLEs this correspondence isn't too surprising, espe- 
cially when one considers that these are the components of the Hamiltonians "perpendicular" 



to the boost that relates the two defining sets of observers. See |T^, |H| for the equivalent 
situation discussed for sets of timelike observers that are boosted relative to each other. 

Less similar are the "boost-dependent" terms which we might expect to be different. 
Starting again with the electromagnetic term, we see that the two are functionally equivalent 
(up to the lapse dependence and factor of two) though their Coulomb potentials are defined 
with respect to different normals. Most different though are the remaining two geometric 
terms which may be considered the key terms in the expressions since they are all that 
remains when the angular momentum and electromagnetic terms vanish. depends on 
the acceleration of the null congruence of generators of A while Hb depends on extrinsic 
curvature of the in S^. The two are not equivalent, for the analogue of the extrinsic 
curvature of Bf is the expansion ^ of A^, which is set to zero by the defining conditions for 
non-expanding horizons. Further, the corresponding acceleration n°'aa on B appears only 
in the stress tensor s"^ which in turn only shows up in calculating the rate of change of the 



CQLE Hamiltonian if Bt is not static []T2|, [14 



Not surprisingly given the difference in the definitions of the two, their numerical values 
differ as well. For simplicity in the calculation I'll drop the electromagnetic field (and the 
associated complications of gauge choice) as well as the angular momentum (since it con- 
tributes in the same way to each term) and for this comparison work in the Schwarzschild 
spacetime. Then, based on the calibration of the previous chapter 

Ha = m, (52) 

where m is the mass parameter for the black holes. Unfortunately, it is not quite so straight- 
forward to give a value for Hb as there are a variety of surfaces B, reference terms e^, and 
lapses N that may be chosen. We will choose the so-called canonical QLE for our compar- 
ison by setting N = 1 (and so normalizing to length —1) and further choose spherical 
Bt surfaces so that we may use the standard embedding reference terms 0. Finally, we will 
choose to calculate the CQLE for a constant radius surface B outside of A and then consider 
the limit as S — A (we could equally well consider a timelike surface B intersecting A at 
At and then take the limit as B becomes null. The result is the same and dealt with in some 



detail in [T^, |Tg). Then 

H^^^^ = r{l- ^l-2m/r), (53) 

which goes to m in the limit r — oo and then monotonically increases to 2m as r — >■ 2m 
(though of course B goes null at 2m and so Hb isn't properly defined right on the horizon). 
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In the literature, this is usually taken to indicate that — m units of gravitational potential 
energy live between the horizon of the black hole and infinity. From this point of view then, 
the fact that the NQLE measures a mass of m at the horizon could be taken as indicating that 
it not only includes the external electromagnetic "hair" but also includes the corresponding 
gravitational "hair" . 

There is also a difference in the reference terms needed to normalize each of the QLEs. 
Namely the reference term on B is allowed to be any functional of 7^/3 and 7f while i^A 
as we have seen is severely restricted in the form that it is allowed to take. This difference 
is perhaps not so surprising when we examine the variational calculations more closely. 
While the mechanics of the calculation are more-or-less identical, there are some significant 
differences in the boundary terms. Recall that in order for the variational principal to be 
well defined for a timelike boundary, the full (induced) metric on that boundary must be 
fixed and unvarying. By contrast, for a null boundary, one need only say that it is a non- 
expanding horizon. The size of cross-sections on any non-expanding horizon is fixed in 
time for a given horizon, but variations are allowed that will move to different horizons with 
different areas, charges, and angular momenta. Because there are fewer fixed characterizing 
terms for an isolated horizon, there is correspondingly less freedom to add constants onto the 
action (and therefore the Hamiltonian), to the extent that dual requirements that there be a 
consistent Hamiltonian evolution and that H/^ vanish when A/^, Ja, and Qa — ^ essentially 
fixes what boundary terms can be added to the Hamiltonian in the null case. By contrast, 
with the entire history and future of the boundary evolution defined in the time-like case, 
there is more freedom to add boundary terms. As long as they are functionals of the boundary 
metric 7q,^ they do not affect the Hamiltonian evolution. Thus, in this case there is an 
indeterminacy in the definition of the QLE that is not found for the NQLE on isolated 
horizons and one is free to choose "reference" terms appropriate to a given situation (for 
example the embedding reference terms of 0] versus the intrinsic counterterms of ||24|| ). 

It is of interest too, to compare the thermodynamics of the two formalisms. To do this 
however, we first need to note that the Brown- York paper [^] makes extensive use of the 
path-integral formalism of gravity to put forward a theory of thermodynamics including 
derivations of the entropy and temperature. By contrast, in this paper we have really been 
looking at isolated horizon mechanics. The difference is that while we have tacitly assumed 
the surface-gravity/temperature and surf ace- area/entropy identities, they cannot be proved 
within the classical framework used (but see for a canonical quantum gravity calculation 
of the entropy of an isolated horizon). Thus, the aims of the two are slightly different. 
Nevertheless with the variational calculations underlying the two essentially the same it is 
of interest to make the comparison, and so we now do that. 

5.2 Thermodynamics 

Again we start with the similarities between the isolated horizon and Brown- York approaches 
to thermodynamics. For the reader who is not familiar with the Brown- York approach it 
is briefly reviewed in appendix 0. Apart from the fact that they are both Hamiltonian 
based the other way in which they depart from textbook black hole thermodynamics is that 
they seek to study finite regions of spacetime rather than entire spacetimes. While in the 
isolated horizon case that region is contained by the isolated horizon itself, in the Brown- 
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York case, the region contained by a timelike boundary some distance from the horizon. 
Usually that boundary is taken to have a static intrinsic metric (though that is not required 
by the formalism as it is in the null case) and it is assumed that the fields inside of B are 
stationary (though the fields outside of B suffer no such restriction). 

However, though they are closely related in these generalities, there are major differ- 
ences. As has been noted and is discussed in appendix 0, Brown- York thermodynamics 
is formulated with path-integral quantum gravity using the Euclidean instanton approxi- 
mation. This instanton is constructed in the usual way by analytically continuing time to 
imaginary values and then periodically identifying it with a carefully chosen period. As a 
result of this process, the black hole horizon is reduced to a conical singularity in a complex 
manifold, and then further reduced to a smooth point no different from any other once an 
appropriate period is chosen. Thus, in a certain sense, the region being studied is a finite 
part of a "complex Euclidean" manifold rather than a region of Lorentzian space containing 
a black hole horizon. 

That said, the entropy of the region being studied corresponds to the action of the 
instanton (which approximates the full density of states function). Its evaluation depends 
crucially in a balance between the acceleration of (as contained in the stress tensor s"^ on 
the horizon) and the chosen period of the imaginary time coordinate. This balance is fixed 
to eliminate conical singularities, but as a side effect gives the black hole entropy as A/ 4. 

The thermodynamics is then extracted by considering variations of the action that cor- 
respond to shifts through the phase space of black holes. The first law takes the form 

6{A/4) = £x [mV^e] - /5^](5[v^(/."JJ + /?(v^s'^V2)5a,,} , (54) 

derived as equation ( |5UD in the appendix. Here the work terms are defined on the boundary 
B rather than the horizon A. The inverse temperatures (3 come from the lapse and shift 
integrated over the periodic time. The first law takes an integral form in this case because 
the zeroth law doesn't hold, which is to say that (3 and VL are not in general constant over 
B and so cannot be pulled out of the integral. Further, in addition to the work terms that 
we found for RRHs, there is an additional work term of the form s°'^5aab (the stress tensor 
contracted with area element distortions). 

6 Discussion 

Hamiltonian investigations of general relativity have long been popular and in particular the 
recent analyses of isolated horizons are based on just such an approach. Though up to now the 
isolated horizon work has been spinor or tetrad based, the current paper has taken a metric- 
based approach. This will be of interest for those who feel more comfortable with metrics 
rather than connections and has also allowed a comparison of the results with the extant work 
on metric-based Hamiltonians, quasilocal energy, and quasilocal thermodynamics. Further 
the calculations of this paper have been based on direct variations of the Hamiltonian which in 
turn was found by a Legendre transform of the action. This approach then serves to provide 
an alternative viewpoint from the symplectic geometry presentations of the existing isolated 
horizon papers, though, as has been pointed out, the two are computationally equivalent. 
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We have seen that the isolated horizons fit very nicely into the metric-based formalism 
with the required quantities that we want to fix naturally occuring as boundary terms in the 
variational calculations. In particular the expansion of the horizon and the electromagnetic 
matter flow across it are naturally fixed, and these of course are key characteristics of non- 
expanding horizons. In agreement with the earlier tetrad and spinor approaches, we saw that 
in order for the variations of the action to vanish across the phase space of RRH-containing- 
solutions to the equations of motion, the time evolution vector field T" must be functionally 
dependent on A^, Ja, and Qa which in turn means that the secondary quantities ki, Q^, 
and must also be dependent on these invariants. As noted, this is a departure from 
ADM-type Hamiltonian treatments of general relativity where the time-flow vector field is 
an independent background structure and does not appear in the four-dimensional covariant 
form of the action. 

We proposed an if^ boundary term defined as a surface integral on At which gives 
rise to a Smarr-type formula for rigidly rotating horizons. Together with the first law of 
RRH mechanics (which arises from the Hamiltonian formalism whether or not we postulate 
that i^A takes this form), this formula then places a fairly strong restriction on the allowed 
functional forms of Ke, Qy,, and The standard K£, Q^^, $£, and energy in Kerr-Newman 
spacetime satisfy these functional forms, and so may be used to calibrate these quantities 
across phase space. It is important to remember however, that even without this calibration 
the Smarr formula and laws of isolated horizon mechanics continue to hold. 

It is perhaps interesting to note that throughout this work, the assumption that (p°' was 
a Killing vector of the RRH was only very weakly used. In fact, the only purpose of the 
rotational symmetry assumption was to select out a vector field that would not be affected 
by the variations, and so could be commuted with 6 in the derivation leading up to the first 
law (equation (|5I1|)). If an alternate method could be found for assigning vector fields in 

across the phase space, then one could easily reformulate all of this in terms of weakly 
isolated horizons rather than RRHs. 

The angular momentum arising from the Hamiltonian has previously been shown to be 
closely linked to the equivalent Komar angular momentum in [|l^. Here, we have also shown 
that our proposed Hamiltonian energy is similarly related to the Komar energy integral. 
In both cases, the purely gravitational parts of the derived quantities exactly match the 
corresponding Komar integral, while the electromagnetic parts do not. In turns out though 
that for stationary spacetimes, those electromagnetic parts serve to make the quasilocal 
quantities defined on the horizon equal to the ADM/Bondi energy defined at infinity. Thus, 
in some sense these terms include the electromagnetic (and gravitational depending on your 
point of view) "hair" of the hole in the quasilocal energy/ angular momentum measured at 
the surface. 

This result is in contrast to the Brown- York quasilocal energy which is based on essentially 
the same calculations as the current one, though with timelike instead of null boundaries 
(in fact the non-orthogonal version of their calculations has been been reproduced on our 
outer boundary B). For their CQLE the energy contained by a surface just outside the 
horizon of a Schwarzschild hole is equal to twice the mass of the hole, while for the NQLE 
defined here it is equal to the mass itself. The CQLE, of course, may also be defined off the 
surface, in which case we find that it monotonically decreases to m at infinity. This is often 
interpreted to mean that is it measuring — m units of gravitational potential energy outside 
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the horizon. Taking that point of view, one can then reconcile the two claims if one thinks of 
the NQLE as including the "gravitational hair" just as it also includes the electromagnetic. 
Then the CQLE could be thought of as measuring the bare mass not including the attendant 
gravitational fields outside the inner boundary. 

Again given the similar mathematical machinery underlying the two approaches it is 
natural to examine the thermodynamics/mechanics to which they give rise. Though both 
the Brown- York and isolated horizon analyses are quasilocal in nature and Hamiltonian 
based, it was seen that in almost all other respects they are quite different. In the first place 
of course. Brown- York aim to provide a theory of thermodynamics while the isolated horizon 
work is less ambitious and just seeks laws that are analogous to the laws of thermodynamics. 
No attempt is made to prove the entropy/area and surface gravity/temperature associations. 
Second, the energy and angular momentum in the Brown- York approach are evaluated on 
a timelike surface apart from the horizon while the corresponding quantities for RRHs are 
calculated directly on A. Third while the zeroth law holds for RRHs it does not for the 
timelike boundaries except in the static, spherically symmetric case (see 0). Fourth, the 
actual use of the variational calculations in the two formalisms is quite different. The first 
law of RRHs essentially arises from a direct variation of the Hamiltonian (identified with 
the energy). By contrast, the Brown- York approach uses the variational principle and path- 
integral gravity in the Euclidean approximation to calculate the total entropy of the region 
of spacetime as the action of an appropriate instanton. Variation of that action the gives 
rise to the first law. 

It would be interesting to investigate more closely what restrictions the enforcement of 
strict Hamiltonian evolution might have on the allowed reference terms of the Brown- York 
action. Up to now they have only been considered the Lagrangian formulation. As has been 
noted however, the situation for the outer boundary terms is not the same as for the horizon 
terms. In particular, on the inner boundary the general-covariance of the action is broken 
by the inclusion of a time-flow dependent boundary term and as a result the time flow is 
restricted by the Hamiltonian formalism. This is not the case at the outer boundary where 
the general covariance has been preserved. Nevertheless, it might be worth thinking on this 
issue more closely to see if any restrictions could be placed on the allowed reference terms. 
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A Brown- York thermodynamics 

In this appendix we summarize the Brown- York approach to thermodynamics (for more 
details see 0. For simplicity we will ignore the electromagnetic field. The quasilocal region 
M is taken to be same as that which we have considered in this paper. That is, a timelike 
outer surface, a null inner surface, and initial and final surfaces Si and S2 that are spacelike 
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everywhere except right at their intersection with A. However, the fohation (and therefore 
Si and S2 as well) is taken to consist of stationary time slices and the time flow vector T" 
is taken to be the stationary Killing vector and an element of TA on that surface. Then 
we must have = on A so that the surface will be null and further that surface will be 
non-expanding because is the stationary Killing vector. Further one assumes that the 
coordinate system is corotating with the horizon, so /N vanishes there. Meanwhile, at the 
outer boundary, T" is assumed to lie in TB, so that boundary is stationary as well. Further 
on 5, microcanonical boundary conditions are imposed. That is e, J"", and are fixed 
and therefore the quasilocal energy, angular momentum, and area of the outer boundary 
are fixed. Thus, looking at the system as a whole, it is isolated from the outside with fixed 
energy, angular momentum, and surface area and contains a non-expanding horizon. 

By analogy with non-gravitational physics, the path integral formulation of quantum 
gravity then says that we can find the density of states of the system by evaluating the path 
integral 

v{ej.a)= JvHpe''-/\ (55) 

where the integral is over all possible spacetime configurations Hp (not just solutions to the 
equations of motion) that are periodic in coordinate time period P = t2 — ti and which 
have a fixed data e, ja, and aajj on the outer boundary. Because of these microcanonical 
boundary conditions, the entropy of the system is approximately equal to log v{e,j, a). 

The catch, of course, is that no one knows how to evaluate such gravitational path 
integrals exactly. Instead, the steepest descent approximation is used. That is, it is assumed 
that the entire integral considered above may be approximated by 

vieJ,cr) = e''*/\ (56) 

where /* is the microcanonical action of a complex and periodic solution to the equations of 
motion that also satisfies the required boundary conditions. Such a solution may be easily 
constructed by analytically continuing —>■ —iN and — > —iV"' for the original solution 
that we are studying. Now after this rotation, A is no longer null. In fact because V^/N —>■ 
there it actually becomes Euclidean. Further, because = the complex spacetime closes 
up at A and a conical singularity appears unless the time period P is carefully chosen. 
Specifically, we must have 

— = n^'VaN = lim 47riV((Tafls"'^). (57) 

Thus, the microcanonical action should be one that fixes the lapse, shift, and pressure s°'^aai3 
on the inner horizon and e, ja, and (Jap on the outer. From the considerations of section ^ 
it is easy to see that the appropriate action is 

Imc = { momentum and constraint terms } (58) 

+ j^^ dt £x^\N{e-^) - V'ja - -js'^'^abj . 

Then, for a stationary complex solution to the Einstein equations that is constrained so that 
N = 0, V"-/N = 0, and equation (|57D holds, this is equal to —iAA/4:. Thus the standard 
entropy result is obtained. 
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The first law is then derived by considering 5Imc- Again leaning on the results from 
section ^ and now assuming a rotational symmetry, we have 

5{A/A) = [ d^x {pS[^/^E] - /31]5[v^<^"jJ + /3(v^s''V2)5(7afe} , (59) 

where we've altered the notation slightly to match that used in this paper. (3 = P is the 
inverse temperature and it comes out of the time-integrated lapse and shift terms. 
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